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Abstract 

Two-dimensional almost-Riemannian structures are generalized Riemannian structures on sur- 
faces for which a local orthonormal frame is given by a Lie bracket generating pair of vector 
fields that can become collinear. Generically, there are three types of points: Riemannian points 
where the two vector fields are linearly independent, Grushin points where the two vector fields 
are collinear but their Lie bracket is not, and tangency points where the two vector fields and 
their Lie bracket are collinear and the missing direction is obtained with one more bracket. 

In this paper we consider the problem of finding normal forms and functional invariants at 
each type of point. We also require that functional invariants are "complete" in the sense that 
they permit to recognize locally isometric structures. 

The problem happens to be equivalent to the one of finding a smooth canonical parameterized 
curve passing through the point and being transversal to the distribution. 

For Riemannian points such that the gradient of the Gaussian curvature K is different from 
zero, we use the level set of K as support of the parameterized curve. For Riemannian points 
such that the gradient of the curvature vanishes (and under additional generic conditions), we 
use a curve which is found by looking for crests and valleys of the curvature. For Grushin points 
we use the set where the vector fields are parallel. 

Tangency points are the most complicated to deal with. The cut locus from the tangency 
point is not a good candidate as canonical parameterized curve since it is known to be non- 
smooth. Thus, we analyse the cut locus from the singular set and we prove that it is not smooth 
cither. A good candidate appears to be a curve which is found by looking for crests and valleys 
of the Gaussian curvature. We prove that the support of such a curve is uniquely determined 
and has a canonical parametrization. 



*This research has been supported by the European Research Council, ERC StG 2009 "GeCoMethods" , contract 
number 239748, by the ANR "GCM", program "Blanc-CSD" project number NT09-504490, and by the DIGITEO 
project "CONGEO". 
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1 Introduction 



A 2-dimensional Almost Riemannian Structure (2-ARS for short) is a rank- varying sub-Riemannian 
structure that can be locally defined by a pair of smooth vector fields on a 2-dimensional manifold, 
satisfying the Hormander condition (see for instance [SJ [TUl [23l [28] ) . These vector fields play the 
role of an orthonormal frame. It can also be defined as an Euclidean bundle of rank 2 on a 2-D 
manifold M and a morphism of vector bundles from E to TM which gives rise to a Lie bracket 
generating distribution. 

Let us denote by A(q) the linear span of the two vector fields at a point q. Where A(g) is 
2-dimensional, the corresponding metric is Riemannian. Where A(q) is 1-dimensional, the corre- 
sponding Riemannian metric is not well-defined. However, thanks to the Hormander condition, one 
can still define the Carnot-Caratheodory distance between two points, which happens to be finite 
and continuous. 

2-ARSs were introduced in the context of hypoelliptic operators (9j[2lJ[22]. They appeared in 
problems of population transfer in quantum systems [15j HU [T7] and have applications to orbital 
transfer in space mechanics [111 [12] . 

Generically, the singular set Z, where A(q) has dimension 1, is a 1-dimensional embedded 
submanifold (see |4j). There are three types of points: Riemannian points, Grushin points where 
A(q) is 1-dimensional and dim(A(<7) + [A,A](q)) = 2 and tangency points where dim(A(g) + 
[A, A](g)) = 1 and the missing direction is obtained with one more bracket. One can easily show 
that at Grushin points, A(q) is transversal to Z. Generically, at tangency points A(q) is tangent 
to Z and tangency points are isolated. 

2-ARSs present very interesting phenomena. For instance, the presence of a singular set permits 
the conjugate locus to be nonempty even if the Gaussian curvature is negative, where it is defined 
(see [1]). Moreover, a Gauss-Bonnet-type formula can be obtained. More precisely, in [U [20] the 
authors studied the generic case without tangency points. In [7] this formula was generalized to the 
case in which tangency points are present. (For generalizations of Gauss-Bonnet formula in related 
contexts, see also [6j [25j [26].) I* 1 [IB] a necessary and sufficient condition for two 2-ARSs on the 
same compact manifold M to be Lipschitz equivalent was given. This equivalence was established 
in terms of graphs associated with the structures. In [19] the heat and the Schrodinger equation 
with the Laplace-Beltrami operator on a 2-ARS were studied. In that paper it was proven that the 
singular set acts as a barrier for the heat flow and for a quantum particle, even though geodesies 
can pass through the singular set without singularities. 

In this paper we consider the problem of finding, at each type of point, a normal form which is 
completely reduced, in the sense that it depends only on the 2-ARS and not on its local represen- 
tation!]] This consists in finding a canonical choice for a local system of coordinates and for a local 
orthonormal frame, i.e., two vector fields F\ and F2 defined in a neighborhood of the origin on M 2 H 
Notice that, the classical normal form in the Riemannian case given by the exponential coordinates 
is not completely reduced in our sense, since it is defined up to a choice of an orthonormal frame 
at the starting point. Rather, for the completely reduced normal forms, which we are looking for 
in this paper, the coordinates and the orthonormal frame should be uniquely determined. 

1 In this paper, an object that depends on the 2-ARS and not on its local representation is called canonical. 

2 To be able to fix completely the system of coordinates and the orthonormal frame, and avoid the problem of 
having quantities defined up to sign, in this paper we assume that the 2-ARS is totally oriented, i.e., both the base 
manifold M and the Euclidean bundle are oriented. See Definition [2] 
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Once a canonical choice of (Fi,^) is provided, the two components of F\ and the two compo- 
nents of F2 are functional invariants of the structure, in the sense that locally isometric structures 
have the same components. Moreover, they are a complete set of invariants because they permits to 
recognize locally isometric structures: if two structures have the same invariants in a neighborhood 
of a point, then they are locally isometric. 

However, one expects that some of these invariants are trivial in the sense that they have always 
the same value. Let us make a rough computation of how many functional invariants we do expect. 
Among the 4 components of the two vector fields, we can fix two by choice of the coordinate 
system (which is defined by a diffeomorphism, i.e., by two functions of two variables) and one by 
choice of the orthonormal frame (which is defined by the choice of an angle, i.e., by a function of 
two variables). Indeed, one of the results of our paper is that under generic conditions, there is 
a canonical way of choosing an orthonormal frame (i 7 !,^) and a system of coordinates wherein 
F\ = (1,0) and F2 = (0, /). Here / is a complete set of invariants (that we call a complete invariant 
since it is just one function). 
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Figure 1: The local representations established in [3J 



Notice that the problem of finding a complete set of invariants is not completely trivial even 
in the simplest case of Riemannian points. See for instance the discussion in [5j [24]. Indeed even 
if one is able to canonically fix a system of coordinates, the Gaussian curvature in that system 
of coordinates is an invariant, but it is not a complete invariant: there are non- locally isometric 
structures having the same curvature (an example is given in Section 13. ip . 

A first step in finding normal forms has been realized in [3], where the local representations 
given in Figure [1] were found. However, the ones corresponding to Riemannian and tangency points 
are not completely reduced. Indeed, there exist changes of coordinates and rotations of the frame 



3 



for which an orthonormal basis has the same expression as in (Fl) (resp. (F3)), but with a different 
function <j) (resp. with different functions tjj and £). 

In order to build the coordinate system to which the local expressions found in |4J apply, the 
following idea was used. Consider a smooth parameterized curve passing through a point q. If the 
curve is assumed to be transversal to the distribution at each point, then the Carnot-Caratheodory 
distance from the curve is shown to be smooth on a neighborhood of q (see [4]). Given a point 
p near q, the first coordinate of p is, by definition, the distance between p and the chosen curve, 
with a suitable choice of sign. The second coordinate of p is the parameter corresponding to the 
point (on the chosen curve) that realizes the distance between p and the curve (see Figure [2]). If 
the parameterized curve used in this construction can be canonically built, then one gets a local 
representation of the form F% = (1,0), F2 = (0,/) which cannot be further reduced. Hence, the 
functions / is a "complete invariant" in the sense above. 




c(-): curve transversal to the distribution 



Figure 2: The construction of coordinates starting from a parameterized curve c(-) :] — e,e[—> M. 
We denote by 7^ the geodesic starting at c(a), parameterized by arclength, entering the region 
where sign = ±1 and such that d(7„ (s), c(] — e, e[)) = s. As the distribution is transversal to c(-) 
at each point, the distance from c(] — e, e[) is smooth. 

For Riemannian points, a canonical parameterized curve transversal to the distribution can be 
easily identified, at least at points where the gradient of the Gaussian curvature is non-zero: one can 
use the level set of the curvature passing from the point, parameterized by arclength (see Section 14.11 
and 14. 1.3[) . For points where the gradient of the curvature vanishes, under additional generic 
conditions, we prove the existence of a smooth parameterized canonical curve passing through the 
point (a crest or a valley of the curvature, see Sections 12.21 14.11 and I4.1.4|) . 

For Grushin points, a canonical curve transversal to the distribution is the set Z. This curve 
has also a natural parameterization, as explained in Section 14. 1\ and was used to get the local 
representation (F2) in Figure [T] (that, as a consequence, cannot be further reduced). 

Concerning the local expression (F3) in Figure HJ in [?| the choice of the smooth parameterized 
curve was arbitrary and not canonical. The main purpose of this paper is to find a canonical one. 
Once this is done, one automatically gets a normal form which cannot be further reduced and the 
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corresponding functional invariant at a tangency point. 

The most natural candidate for such a curve is the cut locus from the tangency point. Never- 
theless, this is not a good choice, as in |13j it was proven that in general the cut locus from the 
point is not smooth but has an asymmetric cusp (see Figure [7]). Another possible candidate is the 
cut locus from the singular set in a neighborhood of the tangency point. The first result of the 
paper concerns the analysis and the description of this locus: in Theorem [T] we prove that the cut 
locus from Z is non-smooth in a neighborhood of a tangency point (see Figure [7] for an example). 
Even if not useful for the construction of the completely reduced normal form, this result is a step 
forward in understanding the geometry of tangency points. 

A third possibility is to look for curves which are crests or valleys of the Gaussian curvature 
and intersect transversally the singular set at a tangency point. The second result of the paper (see 
Theorem [2]) consists in the proof of the existence of such a curve. Moreover, this curve admits a 
canonical regular parameterization. An example of a crest of the curvature at a tangency point is 
shown in Figure O 




-una 



Figure 3: A crest of the curvature at a tangency point for the almost-Riemannian structure on R 2 
having X(x,y) = (1,0), Y(x,y) = (0, y — x 2 ) as orthonormal frame. 

Notice that for 2-ARSs, tangency points are the most difficult to handle due to the fact that 
the asymptotic of the distance to the singular set is different from the two sides of the singular set. 
In [13] the authors gave a description of the geometry of the nilpotent approximation at a tangency 
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point. Also, they provided jets of the exponential map and a description of the cut and conjugate 
loci from a tangency point in the generic case. 

However, tangency points are far to be deeply understood. An open question is the convergence 
or the divergence of the integral of the geodesic curvature on the boundary of a tubular neigh- 
borhood of the singular set, close to a tangency point. This question arose in the proof of the 
Gauss-Bonnet theorem given in [7j. In that paper, thanks to numerical simulations, the authors 
conjecture the divergence of such integral. 

The structure of the paper is the following. In Section [21 we briefly recall the notion of almost- 
Riemannian structure, we define the concept of local representations and describe the procedure 
to build a local representation from a parameterized curve transversal to the distribution. Finally, 
we define the set containing crests and valleys of the curvature. 

Section [3] is devoted to the formal definition of completely reduced normal forms and of complete 
set of invariants as maps that to a germ of a 2-ARS associate a choice of local orthonormal frame 
and a choice of a set of functions which permits to distinguish local isometric structures. 

In Section HJ we state the main results: the construction for each type of point of a canonical 
parameterized curve transversal to the distribution, the description of the cut locus from the singular 
set at a tangency point and the construction of the normal forms and of the invariants. 

The next sections are devoted to the proof of the results. Section [5] is devoted to the description 
of the cut locus from the singular set in a neighborhood of a tangency point. In Section [6j we 
prove the existence of a canonical parameterized curve passing transversally to the distribution 
at a tangency point. In Section [71 we prove Corollary [2] which describes the completely reduced 
normal form and the complete invariant. 

2 Preliminaries 

In this section we recall some basic definitions in the framework of 2-ARS following (4J [7J . 

Let M be a smooth connected surface without boundary. Throughout the paper, unless spec- 
ified, manifolds are smooth (i.e., C°°) and without boundary. Vector fields and differential forms 
are smooth. The set of smooth vector fields on M is denoted by Vec(M). 

Definition 1 A 2-dimensional almost-Riemannian structure (2-ARS) is a triple S = (E,f, (•, •)) 
where E is a vector bundle of rank 2 over M and (•, •) is an Euclidean structure on E, that is, (•, -) q 
is a scalar product on E q smoothly depending on q. Finally f : E — > TM is a morphism of vector 
bundles, i.e., (i) the diagram 

E — ^TM 

7T 

M 

commutes, where ir : TM — > M and tte ■ E —> M denote the canonical projections and (ii) f is 
linear on fibers. Denoting by T(E) the C co (M)-module of smooth sections on E, we assume the 
submodule A = {f o a \ a £ T(E)} to be bracket generating, i.e., Lie q (A) = T q M for every q £ M. 

Denote by Z the singular set of f, i.e., the set of points q of M such that dim(f(£' 9 )) is less than 

2. 
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Definition 2 A 2-ARS is said to be oriented if E is oriented as vector bundle. We say that a 
2-ARS is totally oriented if both E and M are oriented. For a totally oriented 2-ARS, M is split 
into two open sets M + , M~ such that Z = dM + = dM~ , f : E\ M + — > TM + is an orientation- 
preserving isomorphism and f : E\m- ~ > TM~ is an orientation reversing-isomorphism. 

A property (P) defined for 2-ARSs is said to be generic if for every rank-2 vector bundle E over 
M, (P) holds for every f in an open and dense subset of the set of morphisms of vector bundles 
from E to TM, endowed with the C°°-Whitney topology. 

Let S = (E,f, (•,•}) be a 2-ARS on a surface M. We denote by A(g) the linear subspace 
{V(q) | V G A} = f(E q ) C T q M. The Euclidean structure on E induces a symmetric positive- 
definite bilinear form G : A x A — > C°°(M) defined by G(V,W) = (av,crw) where ay,crw are 
the unique sections of E satisfying f o ay = V, f o aw = W. At points q G M where f|# is an 
isomorphism, G is a tensor and the value G(V, W)\ q depends only on (V(q),W(q)). This is no 
longer true at points q where f|p is not injective. 

If ((71,0*2) is an orthonormal frame for (•, •) on an open subset Q of M, an orthonormal frame 
for S on Q is the pair (f o oi, f o 0-2). 

For every q E M and every v £ A(q) define G g (v) = inf{(u, u) q \ u G E q , f(u) = v}. 

An absolutely continuous curve 7 : [0, T] — > M is admissible for S if there exists a measurable 
essentially bounded function [0, T] 9 t 4 u(t) £ Ej(t) such that ^(t) = j(u(t)) for almost every 
t G [0, T\. Given an admissible curve 7 : [0, T] — > M, the length of 7 is 



The Carnot-Caratheodory distance (or sub-Riemannian distance) on M associated with S is defined 

as 



The finiteness and the continuity of d(-, ■) with respect to the topology of M are guaranteed by the 
Lie bracket generating assumption (see [8]). The Carnot-Caratheodory distance endows M with 
the structure of metric space compatible with the topology of M as differentiable manifold. 

Locally, the problem of finding a curve realizing the distance between two fixed points qo,qi G M 
is naturally formulated as the distributional optimal control problem 



where F±,F2 is a local orthonormal frame for the structure. 

A geodesic for S is an admissible curve 7 : [0,T] — > M, such that G 7 ( t )(7(i)) is constant and 
for every sufficiently small interval [ti,t2] C [0, T], 7|[t li t 2 ] is a minimizer of £. A geodesic for which 
G 7 ( t )(7(t)) is (constantly) equal to one is said to be parameterized by arclength. 

If (F\,F2) is an orthonormal frame on an open set O, a curve parameterized by arclength is a 
geodesic if and only if it is the projection on S] of a solution of the Hamiltonian system corresponding 
to the Hamiltonian 




<%o,<7i) 



inf{^(7) I 7(0) = qo,j(T) = 51,7 admissible}. 




H(q,p) = ^{{vFi{q)f + (j>F 2 (q)) 2 ), q G O, p G T*Q. 



(1) 
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lying on the level set H = 1/2. This follows from the Pontryagin Maximum Principle [27] in the 
case of 2-ARS. Its simple form follows from the absence of abnormal extremals in 2-ARS, as a 
consequence of the Hormander condition, see [3]. Notice that H is well-defined on the entire T*M, 
since formula ([T]) does not depend on the choice of the orthonormal frame. When looking for a 
geodesic 7 realizing the distance from a submanifold N (possibly of dimension zero), one should 
add the transversality condition p(0)T 7 ( )iV = 0. 

The cut locus Cutjy from N is the set of points p for which there exists a geodesic realizing the 
distance between N and p losing optimality after p. It is well known (see for instance pQ for a proof 
in the three-dimensional contact case) that, when there are no abnormal extremals, if p E Cutjy 
then one of the following two possibilities occour: i) more than one minimizing geodesic reaches 
p; ii) p belongs to the first conjugate locus from N defined as follows. To simplify the notation, 
assume that all geodesies are defined on [0, 00 [. Define 

C° = {A = (g,p) € T*M I q G N, H(q,p) = 1/2, P T q N = 0} 

and 

exp : C° x [0,oo[-> M 
(A,i) ^vr(e^A) 

where ir is the canonical projection (q, p) — > q and H is the Hamiltonian vector field corresponding 
to H. The first conjugate time from ./V for the geodesic exp(A, •) is 

tconjW = min{i > 0, (A,i) is a critical point of exp}. 

and the first conjugate locus from ./V is {exp(A, t con j(X)) | A G C }. 

2.1 Local representations 

Set Ai = A and A^+i = A& + [A, A&]. Let us introduce the main assumptions under which all the 
results of the paper are proven. 

(HO) (i) Z is an embedded one-dimensional submanifold of M; 

(ii) the points q G M where A.2(q) is one-dimensional are isolated; 

(iii) A 3 (q) = T q M for every q £ M. 

Property (HO) is generic for 2-ARSs (see [4]). 

Definition 3 A local representation of a 2-ARS at a point q G M is a pair of vector fields (X, Y) 
on M. 2 such that there exist: i) a neighborhood U of q in M, a neighborhood V of (0, 0) in M. 2 and 
a diffeomorphism (p : U —> V such that <p(q) = (0,0); ii) a local orthonormal frame (Fi^F?) of A 
around q, such that (p*F\ = X, <~p^F2 = Y , where 99* denotes the push-forward. 

Let us state a result which will be crucial in the following. 

Proposition 1 (|4j) Under the hypothesis (HO), it is always possible to get a local representation 
under the form X = d Xl Y = f(x, y)d y , where f is a smooth function such that one of the following 
conditions holds: /(0,0) 7^ 0, 9^/(0,0) ^ 0, 9^/(0,0) / 0. 
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In the following, we give a procedure which permits to build a local representation of the form 
(d x , f(x, y)dy) starting from a totally oriented 2-ARS and a parameterized curve transversal to 
the distribution. This procedure provides a completely reduced normal form once a canonical 
transversal curve is identified and a proof of Proposition [T^l 

Procedure 1: 

1. Choose any smooth parametrized curve c(-) :] — e,e[—> M such that c(0) = q, 6(a) 7^ for 
q G] — e,e[ and span(c(a)) + A(c(a)) = T C ^M for every a g] — e, e[. 

2. Denote by p :] — e, e[— >■ T*M the smooth map such that, for all a G] — e, e[, 

(a) p(a) G r* a) M, 

(b) p(a)c(a) = 0, 

(c) H(c(a), p(a)) = |, where i/ is defined in formula ([I]), 

(d) p(a)V(a) > (for every a G] — e, e[), if (V(a), 6(a)) is a positively oriented (with respect 
to the orientation of M) pair of vectors applied in c(a). 

Remark that the map p is unique once c(-) is fixed. 

3. Define £ : M? — > M as the map that associates with the pair (x, y) the projection on M 
of the solution at time x of the hamiltonian system on T*M associated with H with initial 
condition (c(y), p(y)). The map £ is a local diffeomorphism from a neighborhood V of G R 2 
to a neighborhood U of g G M and preserves the orientation. Define (p as the inverse of the 
restriction of £ to these neighborhoods. 

4. Let <7 G r(i£|[/) be the unique section of norm one such that </?*(f o cr) = and p such that 
(a, p) G r(£'|t/) 2 is a positively oriented orthonormal frame of E\jj. 

5. Define X and Y in Vec(M 2 ) by X = y?»(f o a) and F = ^(f o p). 

6. Define / : V C M. 2 -> M by F(x, y) = /(», y)9 y . 

Lemma 1 Procedure 1 can be completed. 

Proof: Step 1 is possible thanks to the fact that A(q) has dimension at least one. 

Step 2 is possible by simple linear algebra considerations. This choice of p corresponds to the 
initial condition of a geodesic (solution of the Hamiltonian system defined by H) of the 2-ARS, 
transversal to the curve c(-) at the point c(a), hence minimizing locally the distance to the support 
ofc(-). 

In Step 3, the fact that £ is a local diffeomorphism holds true since f§(0), f|(0) are linearly 
independent vectors, which implies that the Jacobian of £ at is not 0. Indeed §f (0) is the initial 
velocity of a geodesic transversal to the curve c (•) when §|(0) = 6(0). As a consequence ip is a 
local coordinate system in the neighborhood of q. 

3 Proposition[T]is formulated for non totally oriented 2-ARS. One can use Procedure 1 to prove it by fixing arbitrary 
local orientations on the manifold and on the Euclidean bundle. 
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In Step 4 the existence of a is guaranteed by the fact that x h-> £(x,y) are geodesies for all 
y G] — e,e[ which implies that is in the distribution A. Once a is defined, the existence of p 
is a trivial fact. 

In Step 5 the existence of X and Y is guaranteed by the fact that 99 is a local diffeomorphism. 

In Step 6, the fact that Y can be written as f(x,y)d y is the consequence of the Pontryagin 
Maximum Principle. Indeed, since the curves x 1— > E(x,y) are geodesies minimizing the distance 
to the support of c(-) and since these curves have the form x t-t (x,y) in the (x, y)-coordinates, 
then the curves x 1— > (x, y) realize the distance between vertical lines. This implies that any 
solution of the Pontryagin Maximum Principle with initial condition transversal to the vertical axis 
{(x,y) I x = 0} should be transversal to the vertical lines for each x. Hence it should annihilate d y . 
As a consequence of the Pontryagin Maximum Principle it should also annihilate the orthogonal 
to d x for the metric. Hence the orthogonal to the metric is proportional to d y . The fact that / is 
smooth is just a consequence of the smoothness of the metric and of 8. ■ 

Remark 1 Notice that Procedure 1 provides 

• a local coordinate system ip around q such that f(q) = (0, 0) and (p preserves the orientation 



• a 2-ARS on a neighborhood of S M 2 having as positively oriented orthonormal frame 
(d x , f(x, y)d y ). One can easily check that, thanks to (H0) ; / satisfies at least one of the 
following conditions /(0,0) / 0, 03-/(0,0) ^ 0, 0^/(0,0) ^ 0. 

A consequence of Lemma Q] is the following. 

Corollary 1 In the totally oriented case, constructing a local representation of the form X = d x , 
Y = f(x,y)d y is equivalent to choose a parameterized curve transversal to the distribution. 

In [4J the following possible forms for the function / in Proposition Q] were found. 

Proposition 2 (|4j) If a 2-ARS satisfies (HO), then for every point q G M there exist a local 
representation having one of the forms 



where <f>, ip and £ are smooth functions such that (p(0,y) = and ip(0) 7^ 0. 

Remark 2 Notice that in the proposition above we do not take account of orientations. 

Definition 4 Under hypothesis (HO), a point q is said to be a Riemannian point if A(q) = T q M , 
Grushin point if A(q) is one- dimensional and A2(q) = T q M , tangency point if A{q) = A2^q) is 
1-dimensional and A^(q) = T q M . 

Local representations for a Riemannian, Grushin and tangency points are given by (Fl), (F2), (F3) 
respectively. 



ofM; 



(Fl) F 1 {x,y) 
(F2) Fi(x,y) 
(F3) F x {x,y) 



A. 

dx ' 

JL 
dx ' 
d_ 

dx ' 



F 2 (x,y) = e^^, 

F 2 (x,y) = xe« a! >v)$j, 

F 2 (x,y) = (y-x 2 4>(x))e^-§-, 
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Remark 3 The local representations (Fl) and (F3) were obtained arbitrarily choosing a parame- 
terized curve transversal to the distribution. Hence, they are not completely reduced in the sense that 
there are isometric structures having local representations of type (Fl) (resp. (F3)), but with dif- 
ferent functions <p (resp. with different functions tp and!;). In other words, the functions appearing 
in (Fl) and (F3) are not invariants of the structure. 

This topic, together with the proof that the local representation of type (F2) is completely reduced, 
is discussed starting from Section^ 

2.2 Crests and valleys of the curvature 

In this section we find an equation whose solutions contain "crests" and "valleys" of the curvature. 
These loci are used in Section U] to construct a canonical parameterized curve transversal to the 
distribution and to obtain invariants of the structure. 

For simplicity let us assume that the curvature A is a Morse function on M\Z, so that its level 
sets are locally either one-dimensional manifolds (at regular points) or isolated points (at maxima 
or minima) or the union of two one-dimensional manifolds which transversally intersect (at saddle 
points). 

Let us think to A as a function describing the altitude of a mountain. Roughly speaking, crests 
and valleys of A are the loci where the distance among the level sets of A has a local maximum 
along a level set of A, see Figure 0]A. To distinguish among crests and valleys one should consider 
integral curves of —VK. Here VK denotes the almost-Riemannian gradient of A, i.e., the unique 
vector such that G(VK, •) = dK(). Locally, integral curves of —VK diverge from a crest and 
converge to a valley, see Figure EJ3. 

To find an equation satisfied by crests and valleys, let C be a level set of A. For simplicity, 
assume that C is a one-dimensional manifold at q. The point q belongs to a crest or a valley of A if 
||VA~|| 2 := G(VK,VK) has a local minimum on C at q. A necessary condition for this to happen 
is that V||VA|| 2 and VK are collinear at q. Namely, crests and valleys are loci lying in the set 



where (VA) is a non- vanishing vector satisfying GiVK, (VA) ) = 0. Notice that this set contains 
not only crests and valleys but also anti-crests (i.e., loci in which ||VA|| 2 has a local maximum 
along a level set of A and integral curves of —VK diverge from these loci) and anti-valleys (i.e., loci 
in which ||VA|| 2 has a local maximum along a level set of A and integral curves of —VK converge 
to these loci). Moreover the set ([2]) contains critical points of A and ||VA|| 2 . 

For 2-ARSs, A is not defined on the singular set Z, in particular it is not defined at tangency 
points. Since below we look for a canonical curve passing through a tangency point, it is convenient 
to look for it inside the set 4 defined in the following. 

Definition 5 Let S = (E,f, (-, •)) be a 2-ARS. Define 

+ :={qeM\Z\ G(V||VA|| 2 , VK^)\ q = 0} U {q E M \ q is a tangency point ofS}. 

Remark 4 Notice that, for a 2-ARS defined by an orthonormal frame of the form (d x , f(x,y)d y ), 
at Riemannian points the metric and the curvature are given by (see for instance J^) 



{gGM|G(V||VA|| 2 ,VA ± )| (? = 0}, 



(2) 




A = 



2(aj) 2 + fdlf 
P 



(3) 
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crest or valley of K 


crest valley 


K = const^^ 

A 


integral curves of — VK integral curves of — VK 

B 



Figure 4: Crests and valleys of the curvature 



3 Completely reduced normal forms and invariants 

In the following, 2-ARSs are assumed to be totally oriented unless specified. Recall that under the 
generic conditions (HO) we have three types of points: Riemannian (R for short), Grushin (G for 
short) and tangency points (T for short). 

In the following, we need to distinguish two types of tangency points. 

Definition 6 Let q be a tangency point. Let us orient Z as boundary of M~ (see Figure^). We 
say that the tangency point is of type T® (resp. T & ) if the distribution is rotating positively (resp. 
negatively) along Z at q. See Figured 

Remark 5 Notice that the distinction between T® and T® was used also in HE/ , to obtain 
respectively a Gauss-Bonnet Theorem for 2-ARSs and a classification of 2-ARSs w.r.t. Lipschitz 
equivalence (T® corresponds to a contribution r q = — 1 and T® corresponds to a contribution 
T q = 1, where r q is defined in ^). 
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Figure 5: The two different types of tangency points 



In the following we have to make additional generic conditions to treat Riemannian points. 

HA. On M \ Z the Gaussian curvature is a Morse function (i.e., at points where VK = 0, the 
Hessian is non-degenerate). Moreover, if q is such that VK = 0, then the Hessian of K at q 
computed in an orthonormal system of coordinates has two distinct eigenvalues. 
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It is a standard fact that this condition is generic. It ensures that VK vanishes only at isolated 
points and that at these points the curvature has i) a local minimum, ii) a local maximum or iii) 
a saddle. In case iii) the curvature has one crest and one valley intersecting transversallyQ The 
condition on the eigenvalues of the Hessian implies that: in case i) the curvature has one valley and 
one anti-valley intersecting transversally; in case ii) the curvature has one crest and one anti-crest 
intersecting transversally. 

To be able to fix an orientation on crests and valleys we need a higher order condition. 

HB. Assume HA. If c(-) :] — e, e[— > M is a smooth curve parameterized by arclength the sup- 
port of which is contained in a crest or a valley of K and such that Vif(c(0)) = then 
^K(c(t))\ t=0 ^0. 

In the following we call (HI) the collection of HA and HB. We need to treat separately two 
types of Riemannian points. 

Definition 7 Assume (HI) and let q G M\Z. We say that q is a Riemannian point of type Rl 
ifVK(q) / and of type R2 ifVK(q) = 0. 

Recall that a pair of vector fields (X\,X2) defined in a neighborhood of the origin of M? for 
which dim(Lieo{Xi, X2}) = 2 defines a totally oriented 2-ARS, namely the one for which (X\,X2) 
is a positively oriented orthonormal frame and the orientation on M 2 is the canonical one. 

Definition 8 Let Q be the set of germs of totally oriented 2-ARS verifying (HO) and (HI). 

Let O be the set of germs at the origin of pair of vector fields on M 2 such that if o G O then (i) 
o is Lie bracket generating, (ii) the corresponding germ S Q of totally oriented 2-ARS belongs to Q. 

We say that two germs of totally oriented 2-ARS are p-isometric if they are isometric and the 
isometry preserves the orientation of the base manifolds and of the vector bundles^ 

Definition 9 (Completely reduced normal form) A completely reduced normal form ( CRNF, 
for short) for totally oriented 2-ARSs is a map N : Q — > O which associates with a germ S of a 
2-ARS the germ o = N(S) at the origin of a pair of vector fields on M 2 such that 

• S and S are p-isometric; 

• N(Si) = N(S2) if and only if S\ and S2 are p-isometric. 

Definition 10 (Complete set of invariants) Let G°° be the set of germs at the origin of C°° 
functions on M 2 . A complete set of invariants of cardinality k for totally oriented 2-ARSs is a 
map I : Q —> (C°°) fc such that I (Si) = I(S2) if and only if Si and S2 are p-isometric. When the 
complete set of invariants is of cardinality one, it is called a complete invariant. 

Remark 6 Notice that once a CRNF is obtained, one can build a complete set of invariants for 
totally oriented 2-ARSs by constructing it on O, since each element of Q is p-isometric to an 
element of O . One can simply take the non trivial functions among the 4 components of the two 
vector fields of o G O as a complete set of invariants. 

4 Here for convenience we consider that if X7K(q) — then q belongs to the crests, valleys, anticrests, antivalleys 
reaching the point. 

5 A local isometry between the base manifolds preserves the orientation on the vector bundles if the pushforward 
of a positively oriented orthonormal frame for the first structure is a positively oriented orthonormal frame for the 
second structure. 
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3.1 On 2-d Riemannian manifolds the curvature is not a complete invariant 



On 2-d Riemannian manifolds, even when a canonical system of coordinates can be constructed 
(for instance, via a canonical transversal curve and following Procedure 1), the curvature written 
in this system of coordinates is not a complete invariant. Consider for instance on ] - ^[xK the 
two Riemannian metrics, 

( 1 \ ( 1 

l 

(x+l) 2 



Both these metrics are written in a system of coordinates built with Procedure 1, by taking as 
parametrized curve the level set of K passing through the origin and parametrized by arclength. 
From one gets for both metrics that 

9 

K 



However, one easily proves that these two metrics are not isometric, even locally. 



3.2 A CRNF by using a canonical parameterized curve transversal to the dis- 
tribution 

To build a CRNF, we find, for each type of point Rl, R2, G, T, a canonical parameterized curve 
transversal to the distribution and we use Procedure 1 to build a local representation. These curves 
are built in Section [4.11 

Let r be the set of germs of smooth curves taking value on a 2-dimensional manifold. For a 
point q of a 2-d manifold, let us denote T q the subset of germs c(-) G T such that c(0) = q. We 
have to build a map x '■ Q 3 $ > c(-) G T such that: (i) c(-) G r g , where g is the base point of 5; 
(ii) c(-) is transversal to the distribution at g; (iii) x is invariant by p-isometries, that is, if eft is a 
p-isometry between 5 and S' then 0(%(<S)) = x(<->')- 

Once x is built, by using Procedure 1, one can obtain a CRNF N x : Q — > O. By Proposition Q] 
we have the following fact. 

Proposition 3 T/ie image of N x is a subset of O of elements of the form (X±, X2) such that 
X\ = (1, 0) and X2 = (0, /), where f is a germ of a smooth function on M 2 . 

Since x is invariant with respect to p-isometries, N x provides automatically a complete invariant. 

Proposition 4 Let I x : Q -)• C°°(IR 2 ,1R) be the map I X (S) = f where f is the function appearing 
as second component of the second vector field of N X (S). Then I x is a complete invariant. 

The next section is devoted to the construction of the map x- Namely, for each type of point 
we build a parameterized curve transversal to the distribution which is invariant by p-isometries. 
We build the map I x and study its image. 
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the transversal curve is the set Z 



M~ 



distribution 



A: Grushin point (G) 



the transversal curve is a crest of the curvature 




the transversal curve 
is a level set of K 



VK 



C: Riemannian points 
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curvature 
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curvature 




Maxima Minima 

D: Riemannian points (R2) 



Saddles 



Figure 6: Canonical parameterized transversal curves at the different type of points. 



4 Main results 

4.1 Looking for a transversal curve 

In this section we build explicitly the map x which associates with the germ of a totally oriented 
2-ARS at the point q a parameterized curve transversal to the distribution at q. By construction 
this curve is canonical. We treat the different types of points G, T, Rl, R2, separately. 

4.1.1 Grushin points 

Let q be a Grushin point. In a neighborhood of q, Z is transversal to the distribution. It is easy to 
see that for any positively oriented orthonormal frame (Gi, G2), the Lie bracket [G\, G2H.Z modulo 
elements in A does not change, and it is not zero close to q. Hence we fix the transversal curve 
to be the parameterized curve c(-) having the singular set Z as support, such that c(0) = q and 
[Gi, Gyi^c) = c'(a) mod A. Notice that with this choice, the orientation of c(-) is such that on its 
right (w.r.t. the orientation of M) lies the set M + , see Figure E]A. This construction is the same 
as the one made in [1] to build (F2), except that this time we take account of orientations. 

Notice that another possible candidate for a transversal curve could be the cut locus from q. 
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However, this choice would require to prove that the cut locus from q is smooth. 
4.1.2 Tangency points 

The case of tangency points is rather complicated. The first candidate as support of a smooth 
curve is the cut locus from the tangency point. Let us recall a result of [13] where the shape of the 
cut locus at a tangency point has been computed. 

Proposition 5 ( [13J ) Let q £ M be a tangency point of a 2-ARS satisfying (HO) and assume 
there exists a local representation of the type (F3) at q with the property 

^(O)+V(O)0 X £(O,O)^O. 

Then the cut locus from the tangency point accumulates at q as an asymmetric cusp whose branches 
are locally separated by Z. In the coordinate system where the chosen local representation is (F3), 
the cut locus is locally 

{(sign(ai)t 2 , Vl"^ 3 + o(t 3 )) | t > 0} U {(sign(a 2 )t 2 , -\^a^ 3 + o(i 3 )) | t > 0}, 

with en = Ci/(tp'(0) + ip(0)d x ^(0, 0)) 3 , the constants Ci being nonzero and independent on the struc- 
ture. 

In general, due to Proposition the cut locus from q is neither smooth nor transversal to the 
distribution at q. 

Another candidate would be the cut locus from Z in a neighborhood of a tangency point. A 
description of such locus is given by the following theorem. 

Theorem 1 Let q £ M be a tangency point of a 2-ARS satisfying (HO) and assume there exists 
a local representation of the type (F3) at q with the property 

a = ^(0)+V(0)S«e(0,0)^0. 

Then the cut locus from the singular set Z in a neighborhood of q accumulates at q as the union 
of two curves locally separated by Z. One of them is contained in the set {y > x 2 ip(x)}, takes the 
form 

{(-1/2^' (0)t + o(t),t + o(t)) | t > 0}, 

and accumulates at q transversally to the distribution. The other one is contained in the set {y < 
x 2 ip(x)} and takes the form 

{(aujt 2 + o(t 2 ),-t 3 + o(t 3 )) \ t>Q}, 

where lo ^ is a constant depending on the structure. This part of the cut locus accumulates at q 
with tangent direction at q belonging to the distribution. 

The proof of Theorem [1] is given in Section \E\ 

As a consequence of the Theorem [TJ in general the cut locus from Z is neither smooth nor 
transversal to the distribution at q. 

Finally, we look for a smooth curve, transversal to the distribution, the support of which is a 
subset of 4k (see Definition [5|) . More precisely, we get the following result proved in Section [6l 



16 




the set * curves through the tangency point 



Figure 7: The singular set (dotted line), the cut locus from a tangency point (semidashed line), 
the cut locus from the singular set (dashed line), and the set 6 (solid lines) for the ARS with 
orthonormal frame F\ = J|, F2 = (y — x 2 — x 3 )-^. In this case the set 4 is the union of three 
smooth curves. Notice that all these curves but one are tangent to the distribution. 

Theorem 2 Let S be a totally oriented 2-ARS on M satisfying (HO) and let q € M be a tangency 
point. There exist e > and a unique smooth parameterized curve c(-) defined on ] — e, e[ which 
satisfies the following properties: (i) c(0) = q, c(0) ^ A(g); (ii) the support of c(-) is contained in 
4k; (Hi) if (X,Y) is the positively oriented orthonormal frame constructed following Procedure 1, 
there exists A > such that 

. [X, [X,Y]]\ (0ty) = -2d y and [d y ,Y]\ m = X8 y if q is of type T®. 

• [X, [X,Y]]\( 0>y) = 2d y and [d y ,Y]\ {m = -\d y if q is of type T e . 
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Notice that with this choice, c(0) points towards M + for T® and towards M~ for T®, see Figure 
[6)3. It is always in the direction where the cut locus from Z accumulates to the tangency point 
transversally to Z (see Theorem [1]). 

4.1.3 Riemannian points of type 1 

Let q be a point such that VK(q) ^ 0. There exists a unique arclength parameterized curve c(-) 
such that c(0) = q, which is supported in the level set {p G M \ K(p) = K(q)}, and along which 
the gradient of K points on the right of c(-) (with respect to the orientation of M). See Figure [6p. 

4.1.4 Riemannian points of type 2 

Let q be a point such that VK(q) = 0. Thanks to (HI), at maxima and saddles, there exists 
exactly one smooth curve 7 C (-) parameterized by arclength with 7 C (0) = q, whose support is a 
crest, and such that J&-if(7 c (£))|t=o > 0. Similarly, at minima and saddles, there exists exactly 
one smooth curve j v (-) parameterized by arclength with 7„(0) = q, whose support is a valley, and 
such that J^-if(7„(i))| 

t=o > 0. When q is a maximum or a saddle, we set c(-) = 7 C (-)- If 9 is a 
minimum we set c(-) = Tu(-). Notice that the support of c(-) is contained in 4k. See Figure [Hp. 

4.2 Completely reduced normal forms and invariants 

In this section, using the canonical parameterized curves built above and applying Procedure 1, we 
give a CRNF for generic totally oriented 2-ARS which provides a complete invariant, i.e., a map I x 
which associates with a germ of a 2-ARS the germ of a smooth function. Finally we characterize 
the image of I x . Recall that C°° is the set of germs of smooth functions on M 2 . 

Corollary 2 (Main Result) Let N x be the completely reduced normal form for totally oriented 
2-ARS obtained from the map \ with Procedure 1. Then the image of N x has the form (X\,X2) 
where X\ = (1,0), X2 = (0,f(x,y)) with f satisfying the following properties. 

[G] If q is a point of type G then f satisfies 

(Ga) /(0,-) = 0. 
(Gb) d x f(0,-) = l. 

In this case we write f G J(tG C <C°°. 

[T] If q is a point of type T® (resp. T e ) then f satisfies 

(Ta) /(0,0) =0. 
(Tb) d x f{0,0) = 0. 

(Tc) d a 2 /(0, •) = -2 for T® (resp. cg/(0, •) = 2 for T© j. 
(Td) d y f(0, 0) > for T© (resp. d y f(0, 0) < for T® ). 

(Te) for y ^ 0, we have G(V||Vi^|| 2 , VK- 1 -)]^^ = which ensures that the vertical axis is 
included in A. 

In this case we write f £ Jfr® C C°° (resp. f G f>e). 
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[Ri] If q is a point of type Rl then f satisfies 

(Rla) /(O, •) = 1 (resp —1) if at q the manifold M and the oriented 2-ARS have the same 
(resp. opposite) orientation. 

(Rib) The second compontent of VK vanishes along the vertical axes. 

(Rlc) The first component of X7K is positive along the vertical axes. 

In this case we write f G C C°°. Writing f = ±e^ x ' y ^ these conditions read 

(Rid) <f)(0, •) = 0. 

(Rle) -2d%<j>(0, y)d x d y <f>{0,y) + 8^(0, y) = 0, for all y. 
(Rlf) 8g0(O, y) - 2d x <f>(p,y)ti%<l>(p,y) > 0, for all y. 

[R2] If q is a point of type R2 then f satisfies 

(R2a) /(O, •) = 1 (resp —1 ) if at q the manifold M and the oriented 2-ARS have the same 
(resp. opposite) orientation. 

(R2b) For y 7^ 0, we have G(X7\\X7 K\\ 2 , VK^-)\( y\ = which ensures that the vertical axis 
is included in <k. 

(R2c) If q is a local maximum for K, then > dyK(0, 0) > d^K(0, 0) which ensures that 
the vertical axis is a crest (and the horizontal one an anticrest). 

(R2d) If q is a local minimum for K, then < dyK(0, 0) < d^K(0, 0) which ensures that 
the vertical axis is a valley (and the horizontal one an antivalley). 

(R2e) If q is a saddle for K, then d%K(0,0) > > d%K(0,0) which ensures that the vertical 
axis is a crest (and the horizontal one a valley). 

(R2f) dyK(0, 0) > which fixes the orientation of the vertical axis. 
In this case we write f G «^r 2 C C°°. 

Remark 7 Notice that the sets Jt»G> $*®> 1*e> Aru Xr 2 are disjoint. 

Moreover their union is not <C°° . This is a consequence of the fact that we have generic condi- 
tions and that f is constructed using a canonical transversal curve. When the function f is obtained 
applying Procedure 1 to any transversal curve, it does not satisfy the conditions given in Corollary 
IE in general. 

5 Proof of Theorem [D the cut locus Cut^ from the singular set 

Z 

In this section, we prove Theorem [1] starting from the local representation (F3). Notice that by 
applying the coordinates change 

y 

x = x, y 



we may assume that "0(0) = 1- For sake of readability, in the following we rename x, y by x,y. 
Since ip(0) > 0, the singular set Z is locally contained in the upper half plane {(x,y) \ y > 0}. 
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Locally, the singular set separates M in two domains {(x,y) \ y — x 2 ip(x) > 0} and {(x,y) \ 
y — x 2 il)(x) < 0}. First, notice that Cut^ D Z = 0, since we are computing the cut locus from 
Z. Second, thanks to hypothesis (HO), the only points of Z where Cut ,2 may accumulate are 
the tangency points, since all other points of Z are Grushin points, where A is transversal to Z. 
Hence, close to a tangency point, Cut 2 is the union of two parts, Cut^ lying in the upper domain 
{(x,y) I y — x 2 tp(x) > 0} and Cut^ in the lower one. 

Applying the Pontryagin Maximum Principle, geodesies for the ARS are projections on M 2 of 
solutions of the Hamiltonian system associated with the function 



H 



1 



that is, solutions of the system 



x 

y 

Px 

Py 



x 2 i){x))e^ x ^f 



Px 

p y ({y 

p 2 (y - x 2 ip(x))(2xip(x) + x 2 i)'(x) — (y — x 2 ip(x))^(x,y))e 2 ^ x ^ 
-P 2 y {y - x 2 ^{x)){\ + {y- x 2 iP(x))^(x, y))e 2 ^. 



(4) 



In addition, a geodesic starting from Z with x(0) = a / 0, realizing the distance from Z and 
parameterized by arclength, must satisfy the transversality condition 

1 



Px (0)=±l, p y (0) = T: 



2aip(a) + a 2 ifj'(a) ' 

As we shall see, the two components of Cut^ have different natures. In the upper domain the 
geodesic starting at a point (a, a 2 ^(a)) and transversal to Z reaches its cut point at a time of order 
1 in \a\. Whereas in the lower domain the geodesic starting at the same point reaches its cut point 
at a time of order 1 in v/|a[. 

Introducing the new time variable s = - where rj > is a parameter, system flU) becomes 



tlx 
ds 
dy 
ds 
dp a 

ds 
dp v 
ds 



x 2 ^{x))e^ x ^f 



-- VPx 

- VPy((y 

'- VPl(y - x 2 if)(x))(2xilj(x) + x 2 %l)'(x) - (y — x 2 4)(x))^(x,y))e 2 ^ x ^ 

3* = -^(2/-xV(x))(l + (2/-xV(x))|(x l2 /))e 2 ^. 



(5) 



In order to make the coordinates x, y, p x and p y dependent on the time variable and on the initial 
condition a, in the following we write them as functions of t and a. In Section [5.11 the parameter 
p x (0) is assumed to be — sign(a), since this implies that the geodesic enters the upper domain. In 
Section T5.2I it is assumed to be sign(a) since this implies that the geodesic enters the lower domain. 

In the following, since we are studying Cut_g close to the tangency point (0, 0), a can be assumed 
as small as we want. 



5.1 The upper part of the cut locus 

We consider the geodesic starting from a point of Z, realizing for small time the distance from Z 
and entering the upper domain, that is, with the initial conditions 

x(t = 0, a) = a, p x (t = 0, a) = — sign(a), 



y(t = 0, a) = a 2 t/j(a), p y (t = 0, a) 



sign(oj 



(6) 



2a*0(a) +a 2 ip'(a)' 
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Computation of jets. For a > 0, choosing rj = a in ([5]) and writing x, y, p x and p y as functions of 
a and s, one can check that if x, y,p x ,p y have orders 1, 2, 0, —1 in a respectively, then the dynamics 
has the same or higher orders. As a consequence, since the initial conditions respect these orders, 
we can compute jets with respect to a of the solution of system ([5]) under the form 

x(s,a) = ax (s) + a 2 xi(s) + a 3 x(s,a) p x (s,a), = p xQ (s) + ap xl (s) + a 2 p x (s,a), 

y{s,a) = a 2 y (s) + a 3 yi(s) + a 4 y(s,a) P y (s,a), = a' l p yQ {s) + p yi (s) + ap y {s,a), 

where x,y,p x ,p y are smooth functions. Using ©, the initial conditions are given by 

a? (0) = 1, xi (0) = 0, p x0 (0) = -1, Pxl (0) = 0, 

y (0) = l, yi (0)=^'(0), Py (0) = h Pyi(0) =-§</>' (0), 

and from system ([5]) we easily get 

x (s) = l-s, xi(s) = 0, yo 0) = 1, yiO) = i//(0), 

whence 

x(t, a) = a — t + a 3 x(a,t/a), y(t, a) = a 2 + a 3 if)'(0) + a y(a, t/a). 
Similarly, for a < one gets 

x(t, a) = a + t + a 3 x(a,t/a), y(t, a) = a 2 + a 3 ip'(0) + a A y(a, t/a). 

Lemma 2 For a small enough and t such that |^| < 2 one gets that 

dx{t,a) 1 dy(t,a) dy(t,a) 

— a >o> — a <0ifa<0, — > if a > 0. 

oa 2 oa oa 

The proof is a direct consequence of the computation of the jets of the geodesies. 

Lemma 3 A geodesic with the parameter a small enough intersects the geodesic with initial con- 
dition a = —a — a 2 ip'(0) + o(a 2 ) at time ti nt (a) = |a|(l + ^aifi'fo) + o(a)). The pair (a,tj n j(a)) is 
unique among the pairs (b,r) realizing the intersection and satisfying ab < 0, r > 0, < |-| < 2 
and < |^| < 2. 

A direct consequence of this lemma is that the cut time of a geodesic starting from Z with initial 
condition a is bounded from above by ti n t{a). 

Proof: Assume a > 0, the proof being the same for a < 0. Let a < and t > be such that 
< I < 2 and < ||| < 2. If y(t,a) = y(t,a) then 

a 2 + aV (0) + o(a 3 ) = a 2 + a 3 ip'(0) + o(a 3 ), 

whence a = —a — a 2 -0'(O) + o(a 2 ). Moreover, a) = x(t, a) implies t = a + ^a 2 ip'(0) + o(a 2 ) and 
the intersection point is 

Xint(a) = -^p-a 2 + o(a 2 ), y int (a) = a 2 + o(a 2 ). (7) 

This, together with the fact that dx< ^^ > 5, proves the uniqueness of the pair (a,ti nt (a)). 
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As far as the existence is concerned, it is not hard to compute that the two fronts corresponding 
to positive and negative parameters transversally intersect at the point (xi nt ,yi n t)- Hence, the jets 
computed above are sufficient to show that the geodesic corresponding to a intersects a unique 
geodesic corresponding to an initial condition of the form a = —a — a 2 ip'(0) + o{a 2 ) at a time of the 
form ti n t(a) = |a|(l + \a^' (0) + o(a)). 

■ 

Lemma 4 The conjugate time of a geodesic is strictly bigger than its cut time. 

Proof: Thanks to the previous computations, the absolute value of the Jacobian of the map 
(t, a) i y (x(t, a),y(t, a)) is 2a + 3a 2 ijj'(0) + a 3 H(a, |) where S is a smooth function. This allows to 
conclude that for |-| < 2 and a small enough, the Jacobian is nonzero, hence t is not a conjugate 
time. ■ 

End of the proof. Let us show that if a is small enough, then the geodesic corresponding to 
a is optimal on [0,ti n t(a)]. By contradiction, assume the geodesic looses optimality at a time 
< t < ti n t(a). 

1. t is not a conjugate time thanks to Lemma [H Hence there exists 6^0 such that the geodesies 
corresponding to a and b intersect at time i and they are both optimal on [0,i]. 

2. i cannot satisfy ||| > 2 since this would imply that the geodesic corresponding to b is not 
optimal until t thanks to Lemma El 

3. If ||| < 2, by the uniqueness property given in Lemma [3l the parameters b and a cannot have 
opposite signs. Nevertheless they can neither have the same sign as Lemma [2] implies that 
two geodesies corresponding to different parameters of the same sign cannot intersect at time 
t. 

Hence t»nt(a) is the cut time of the geodesic corresponding to the parameter a and the cut point is 
given in ([7]). By the transversality of the two fronts (corresponding to initial conditions a and a) 
at time tint, the set Cut^ is locally a 1-dimensional manifold. Moreover, formula (JT]) implies that 
Cutij is transversal to the distribution at (0,0), its tangent vector at (0,0) being (— ip'(0)/2, 1). 

5.2 The lower part of the cut locus 

Reasoning as in section 15-H we consider the geodesic starting from Z, realizing the distance from 
Z and entering the lower domain, that is with the initial conditions 



x(t = 0, a) = a, p x (t = 0, a) = sign(a 

y(t = 0, a) = a 2 ip(a), p y (t = 0,a) = 



2,,,.,, ., _ „ „, _ sign(a) (8) 

2aV>(a) +a 2 ip'{a)' 

Computation of the jets of the geodesies. For a > 0, setting r] = \/a and s = ^, one can check 
that if x, y,p x ,Py have orders in rj higher or equal to 1, 3, 0, —2, respectively, then the dynamics has 
the same or higher orders. As a consequence, since the initial condition respects these orders, one 
can compute jets with respect to r] of the solution of system ([5]) under the form 

x(s,7]) = rjx (s) +rj 2 xi(s) +rj 3 x(s,rj), p x (s,rj) = p x0 (s) + r)p xl (s) + rfp x (s,n), 

y(s,rj) = ri 3 y (s) +rj 4 yi(s) +rj 5 y(s,ri), P y (s,v) = T 2 Py ( s ) + ^Pyii 8 ) + Py(s, v), 
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where x,y,p x ,p y are smooth functions. From the initial conditions ©, we deduce 

z (0) = 0, xi(0) = 1, Px0 (0) = 1, p xl (0) = 0, 
y (0) = 0, y 1 (0) = l, p yo (0) = -|, p yi (0) = 0, 

and using system ([5]), the functions xq, xi, yo, yi,Pxo,Px\,Pyo,Py\ satisfy 



x 

yo 



Pxo 
2 4 

7 Pyo x o 



Pxo = -27 2 Pyo x O 
I Pj/ = 



Xl = Px\ 



Pxi = 7 Py x o(-^P yi Xo + 2py Q y - 6p yQ x xi - 5ap yQ a 



(9) 



I Pvi 



2 2 
7 Pyo x O 



where 7 = e^°'°) and a = t/)'(Q) + f§(0, 0). Thus p yo = — \ and one can prove (see [31 [13]) that 



yo 0) 



2 



3^7 



cn (K + yfis), 



where /C is the complete elliptic integral of the first kind of modulus , and cn , sn and dn 

denote the classical Jacobi functions of modulus -4=. Recall that the Jacobi functions cn,sn are 
4/C-periodic, when dn is 2/C-periodic. 

Denote by xio, yio,Pxio,Py\o the solution of the second system in with a = 0. Define 
91,92,93,94 by 

xi = x w + agi, p xl = Pxio + ag^, 

y\ = Vw + ag2, Pyi = P yi0 + ag 4 . 

It is easy to see that the gi satisfy 



f 54 = 

91 = 93 

92 = 

I 93 = -i7^o(%i + 5x'o), 



-l 2 xU 2 9i +x 2 q ) 



(10) 



and the initial conditions are <7i(0) = 52(0) = 53(0) = 0. Notice moreover that, if 

(X0, yo,Px ,Py , X W ,yi0,Pxl ,Py 10 , 91,92,93) 

is the solution of Q, ([TO]) with initial condition (0, 0, 1, —1/2, 1, 1, 0, 0, 0, 0, 0) then the solution of 
@, CED with initial condition (0,0,-1,-1/2,-1,1,0,0,0,0,0) is 

(-^0, V0, -Px0,Py Q , -X W ,y W , -PxlO,Py 10 ,9l, -92,93), 

which corresponds to the geodesic starting from Z with the initial condition —a < 0. 



Lemma 5 If 5 > and rj ^ are small enough and < |^| < ^= + 5 then — ^ 



dx(t,a) 



> 0. 
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Proof: Assume a > (the computation being the same for a < 0). Then 

" - x (-] - -x (-] + rj (2xi ( - J - -xi (-] J + r/ 2 x r f 77, * 



~u 1 1 ~u 1 I'll ± 1 1 — j. 1 ' I — ' I " 

or] \r] J rj \rj J \ \rj J rj \V / / V V , 

where x r is a smooth function. Now, the function / : u 1— > xq(u) — uxq(u) is such that /(0) = and 



f'(u) = —uxq(u) = -uj 2 Xq(u) > for u G 



0,^ 



Hence, for e small enough, there exists (5 > small enough such that f(u) > e for u &]5, ^= + 5[. 

Therefore, if 5 < l - < ^ + 5, then > e + 77 (2x1 (j) - ^Xi + r/ 2 x r (77, ±) > for 77 

small enough. 

For < I < 5 (possibly reducing (5 and 77), since 2xi(|) — ^xi(^) = 2 for t = 0, we have that 
dz(f,a) - n ■ 

Lemma 6 y4 geodesic with a > intersects a geodesic with a < at time ti nt (a) where 

a = -a+ ^g 2 (^)a^a + o(ay/a), 

tint (a) = ^(Va-Vl ^-je ^a + o(a)). 

Proof: In order to find the expressions given in the statement, we proceed as follows: we fix 
a time to = ^fVo an d we nn d two parameters a > and a < such that the corresponding 
geodesies intersect at time to- Indeed, to is a natural candidate to approximate the intersection of 
the geodesies with initial conditions t}q and — t]q since ^= is the half period of xq. 
We look for a and a by setting 

r /+ = ^ = r /0 + c + r /0 2 + O (r /0 2 ), 

v- = VW\ =m + c-vo + °(vo), 

where c+ and c_ are constants to be found. Remark that this is equivalent to choose to = ^=(77+ ~ 

c+rj+ + o{rj\)) and ry_ = 77+ + (c_ — c + )r]\ + o{r]\). The corresponding geodesic parameterized by 
s are 

x +(s) = r] + x (s) +rf + (x 10 (s) + ag^s)) + o(r/ 2 ), 
V+( s ) = V%yo(s)+r]i(yw(s) + (xg2(s)) + o(r)l), 

and 

x-(s) = -r)-xo(s) + T] 2 (-xio(s) + agi(s)) + o(r/ 2 ), 
j/-(s) = 7£yo(s)+7/i(yio(s)-a52(s)) + o(»7i)- 

Let us estimate the geodesic corresponding to 77+ (resp. rj-) at s+ = (resp. s_ = ^-). One 
computes easily that 

2/C , .. 

s+ = — (1 -c + ?7o + o(r?o)), 

2/C 

s- = —(1- 0-770 + 0(770)), 
V7 
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and 



x + (s + ) 



2 / 2K (2K\ f2IC\\ , 2, 

3 4/C : 4 ( 4tCc+ . ( 2/C \ /2/C 



y + (s + ) = -^ + Vo[-^ T + yio[^=)+a 92 ^—)) + o( Vo ), 

x-(a-) = % (^-c_— -x 10 (^— J +a 5l JJ +0(7/0), 
3 4/C ; 4 / 4/Cc_ x (2K\ (2K 



M *- ) = -"»vf + "»l~ +!,10 Ui" a£ ' 2 U ll+ ° ( " o) 

Hence, these two geodesies intersect at time to for 

ag 2 (^)-2x w (^) 



c+ = ^7 



The intersection point is 



4/C 

a 92 (^) + 2x 10 (^) 
^ 4/C ' 



^nt(to) = %OL ^ ^ + 0(77^), 

I/in*(*o) = ~^^+o(Tjg), 

and the intersection time satisfies 

tn 2/C , ,, 2/C . ._ 

— = —(1-0+770 + 0(7/0)) <— + <5, (12) 

for a small enough and 5 as in Lemma [SJ 

One can compute that the two fronts (corresponding to positive and negative initial conditions) 
are transversal at time ti n t(a) at the point {xintiVint)- Hence the computation with the jets allows 
to conclude that the two geodesies intersect. 

Notice that numerical computations show that 2(ft(^=) + g 2 (-^=) ^ 0. Moreover, thanks to the 
assumption in Theorem [TJ ■ 

Lemma [6] implies that a geodesic with initial condition a looses optimality at t < ti n t(a) which 
is less than ^= + 5 for a small enough. 



Lemma 7 The conjugate time of a geodesic is strictly bigger than its cut time. 

Proof: Thanks to the computations on the jets, one can show that the absolute value of the 
Jacobian of (s,7?) h-» (x(s,rj),y(s,r])) is r] 3 (J (s) + T]Ji(s) + rf J 2 (s, 7/)), where J (s) = x (s)y (s) - 
3yo(s)xo(s), Ji(0) = — 4sign(xo(0)), and J 2 is a smooth function. It was proven in [T3] that Jo 
never vanishes between and s with s > 2JC/^fy. Moreover Ji(0) has the same sign as the function 
Jo on the interval ]0, s[. This allows to conclude that for S > small enough, the Jacobian never 
vanishes on the interval [0, ^= + 5[ which implies that if | < ^7= + 5 and a is small enough then t 
is not a conjugate time for the geodesic with initial condition a. ■ 
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Variations of x(-,a) until ti n t(a). Let us consider a geodesic with a > 0. Since xq satisfies 
^<0on ]0, ^gi(O) = 1, ^(^=) = -1, one can prove that x satisfies §f (s,r/) > on the 
interval [0, ^= + e] and §§(s, ??) < on the interval [^f= + e, ^= + 5] for rj small enough where e is 
a small parameter of order 1 in rj. In particular between t = and t = imt(o), £ is first increasing 
and after decreasing until ti n t(a). 

For a < 0, one can prove the same way that between t = and t = ti nt (a), x is first decreasing 
and after increasing until ti n t(a). 

Estimation of x(-,a) after the first intersection with the x-axis. Consider a geodesic with 
a > 0. Call s a the first time s such that y(a, s) = 0. For any A > 0, one can compute that 
y(s, Tj) = T) + o(?? 4 ) > for 77 small enough and any s £ [0, A77] . Hence s a > \i] for i] small enough. 

Since x(Xrj,rj) = rf(l + A) + o(r] 2 ), fixing A > a and thanks to the previous 

considerations on the variations of x, the minimum of the ^-coordinate for t G [Tjs a ,ti n t(a)] is 
attained at t = ti n t{a). 

For a < 0, one can prove the same way that the maximum of the x-coordinate for t £ [rjs a , imt(a)] 
is attained at t = ti n t(a). 

Notice that for any geodesic of the 2-ARS the y-coordinate is monotone thanks to the normal 
form (F3). Indeed, one easily proves that horizontal lines parameterized by arclength are optimal 
geodesies. 

End of the proof. Now, we have all the ingredients to conclude. 

1. A geodesic cannot loose optimality by reaching its conjugate locus thanks to Lemma [71 

2. Assume that two geodesies with initial conditions a 7^ a 1 of the same sign loose optimality by 
intersecting one each other. They should intersect before U n t(a) and tj n t(a'), i.e., for t such 
that < |^| < ^7= + 5 and < |^-| < ^= + 5. Hence Lemma [5] applies, which leads to a 
contradiction with x(t,a) = x(t,a'). 

As a consequence a geodesic with parameter a looses optimality by intersecting another 
geodesic with parameter a' such that aa' < 0. 

3. Thanks to the monotonicity of the y-coordinate, two geodesies with initial conditions a and 
a 1 of opposite sign can intersect only in the half plane y < 0. 

In the following, we assume 

tha t a ^SJ^Oi < 0, i.e., x (ti n t (a), a) < for a small 
enough, the proof being the same in the opposite case. 

4. Let us consider the geodesies corresponding to a > and a < as in Lemma [6J For a 1 such 
that < a' < a we have that 

x(t, a) > x(t int (a), a) > x(t int (a), a) = x(t int (a), a) > x(t, a) 

for t < ti n t{a!) such that y(t,a) < 0, where the first and the last inequalities follow from 
the estimations of x given above. This implies that the geodesic corresponding to a' cannot 
intersect the one corresponding to a before loosing optimality. 

5. Thanks to Lemma O if a' > a then x(t, a) < x(t, a') for t < tj n t(a) = ti n t(a). Moreover, the 
estimations of x given above imply that x(t, a) < x(t,a) for t < tint (a) such that y(t,a) < 
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and y(t, a) < 0. Hence the geodesic corresponding to a' cannot intersect the geodesic 
corresponding to a before U n t(a). 

Finally we conclude that any geodesic with parameter a < looses optimality by intersecting 
the geodesic with parameter a > 0, which implies that the geodesic with parameter a looses 
optimality at the same time. 

Since the two fronts are transversal at the intersection of the two geodesies corresponding to 
the initial conditions a and a, the lower part of the cut locus is locally a 1-dimensional manifold. 
Together with the formulae (jlip . this implies that the lower part of the cut is half a cusp tangent 
to y = 0. The set Cut^ for an example of 2-ARS is portraited in Figure El 



y 


Z 




^-—^^ X 








geodesies 


wave fror^^ 

cut locus 





Figure 8: The set Cut 2 for the ARS with orthonormal frame F\ = J|, F2 = (y — x 2 — x 3 )-^ 



Remark 8 In U3]j the description of the cut locus to a tangency point was given. In that paper, 
only the existence of the intersection point is actually proved. The same arguments as in Section^ 
can be used to prove that the intersection really corresponds to the cut locus. 



6 Proof of Theorem 



Let us start by proving the existence of a subset of 4ft passing through the tangency point q and 
satisfying the following conditions: i) it is the support of a smooth curve, ii) it has a tangent 
direction which is transversal to the distribution at q. 

Choose a local representation of the type (F3). By construction, K is well defined outside the 
singular set Z. The set 4k \ Z is implicitly defined by the equation 



G(V(||VK|| 2 ),(V^) = 0. 



(13) 



Computing the left hand side of equation (fT3]) , using the expression of the curvature ([3]) , we find 
that 

|2 , , n „ u , e 2 ^ x ^h(x,y) 



G(V(\\VK\\'),(VK)^ 



(y — x 2 ip(x)) 8 ' 
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where h is a smooth function. Hence, equation (|13p is equivalent to h(x, y) = 0. The development 
of h at the point (0, 0) is 

6 

h(x, y) = a;(y 4 (10V(0) 2 x + y(3^'(0) - 20^(0, 0)^(0))) + £ a,(x, y)x i y & - i ) , 

i=0 

where w is a nonzero constant and ai are smooth functions. Let us show that there exists a smooth 
function b : I — > R defined on a neighborhood / of such that after the coordinate change 

x = 10V(0) 2 x + 2/(3^(0) - 20^(0, 0)V(0)) - b(y)y 2 , y = y, 

we have h(x(x,y),y(x,y)) = x/i(x,y). In the new coordinate system, we have 



h(x(x, y),y{x, y)) = u(y A x + F(x, y)), where F(x, y) = b(y)y 6 + ^ di(x(x, y),y)(x(x, y)fy e 



i=0 

In order x to be factorizable in F, we require that F(0, y) = 0. Since F(0, y) = y 6 R(b(y),y), where 



R(b(y),y) = b(y) + £ ^^^ (-3^(0) + 20,£(O, 0)^(0) - b{y)y)\ 

it follows that F(0, y) = if and only if there exists a smooth function 6 defined on a neighborhood 
of such that R(b(y),y) = 0. Let b = ~ELo (-3^(0) + 20^(0,0)^(0))*. Then, since 

(6, y) I—)- y) is smooth, R(b, 0) = 0, dbR(b, 0) = 1, by the implicit function theorem there exists 
a smooth function b(y) with the properties above. Therefore, coming back to the (x, y) coordinates 
we have shown that 

h(x, y) = w(10z + 2/(3^(0) - 20,£(O, 0)) + b(y)y 2 )(y 4 + F(x, y)), 

where F is smooth function of order 5 in (x, y) and 6 is the function built above. The last equation 
implies that the set C = {(x, y) \ Wx+y(3?p' (0)— 2d x ^(0, 0))+b(y)y 2 = 0} is a connected component 
of the set 4*. it passes through (0,0), it is smooth at (0,0) and its tangent line at (0,0) is 

-(25,e(0,0)V(0)-3V'(0))y, 



iov>(o) : 

that is transversal to the distribution at (0,0). 

Moreover, since F has order 5 in (x, y), any curve contained in the set {(x, y) | h(x, y) = 0} but 
not in C must have a tangent line at (0, 0) belonging to the distribution at (0, 0). 

Requiring that the tangent direction to C at (0, 0) is vertical we get the condition 

20^(0,0)^(0) -3^(0) =0. (14) 

As a consequence, any curve transversal to Z at q contained in 6 should have C as support. 

Now, let us prove that there is a canonical parametrization of C. Let us choose any parameteriza- 
tion of C and construct (x, y) and (X, Y) as in Procedure 1. Then, since dim(A(g)) = dim A2(q) = 1 
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and dimA 3 (<2) = 2, one gets that F(0,0) and [X,Y](0,0) are horizontal and [X, [X,Y]](0,0) is not 
horizontal, which is equivalent to /(0,0) = f x (0,0) = and /^(0, 0) ^ 0. 

Assume ./^(O, 0) > 0. Denote by c(-) the parameterized curve whose support is C, such that 
c(0) = q, and satisfying c(y) = l/2f xx (0,y)d y in the coordinates. With this choice for the 

parameterization, if (x, y) and (X, Y) denotes the coordinates and the local representation defined 
by c(-) via Procedure 1, then [X, [X, Y]](0, y) = 2dy. Equivalently, if Y = f(x,y)dy we have 
f xx (0,y) = 2. An easy computation shows that choosing the parameterization y \— > c(—y), the 
condition f xx (0,y) = 2 is still fulfilled. 

If fxx(0, 0) < 0, the same arguments prove that, denoting by c(-) the curve whose support 
is C, such that c(0) = q and satisfying c(y) = — l/2f xx (0,y)d y in the (x,y) coordinates, then 
/xx(0, y) = —2. Moreover, the last condition does not change when choosing the opposite orientation 
for c(-). 

The curve c(-) defined above satisfies (i), (ii) and, depending on the sign of /^(O, 0), the first part 
of (iii). It is uniquely determined up to orientation. One can prove easily that if fydy = [dy,Y] = 
Xdy then, changing the parameterization of c(-) for y h-> c(—y) one gets fydy = [dy, Y] = —Xdy. 

In Theorem [2l we fix the orientation of c(-) (and hence we fix c(-)) by asking that when 
jW(0,0) > then jg(0,0) < and when /^(0,0) < then / 5 (0,0) > 0. It is not hard to 
prove that the first situation corresponds to a point of type T e and the second to a point of type 

7 Proof of Corollary [2] 

Grushin points. Equation (Ga) is a consequence of the fact that the vertical axis x = is 
contained in Z and hence the distribution has dimension one for x = 0. Equation (Gb) expresses 
the fact that [X±, X2] = d y along the vertical axis, by construction. 

Tangency points. Denote by (x, y) the coordinate system and by (X, Y) the orthonormal frame 
given by Procedure 1 when the curve c(-) is the one given by Theorem [2j The fact that A(g) and 
A2(g) have dimension one implies (Ta) and (Tb). Equations (Tc) and (Td) are implied by (iii) of 
Theorem [2l Equation (Te) is a direct consequence of (ii) of Theorem [2l 

Riemannian points of type 1. Properties (Rla), (Rib) and (Rlc) are direct consequences of 
the discussion given in Section 14.1.31 

Equation (Rid) comes from property (Rla). For what concerns (Rle) and (Rlf), let us compute 
VK using the equality ([3]) 

VK(x,y) = {d 3 x ^(x,y)-2d x ^(x,y)d 2 x ^x,y),e 2 ^ (-2d x <i>(x, y)d x dy<f>{x, y) + d^dy^y)) ). 

Applying (Rib) and (Rlc) to this equality, one gets (Rle) and (Rlf). 

Riemannian points of type 2. The properties (R2a) and (R2b) are direct consequences of the 
discussion given in Section 14.1.41 Inequality (R2c) follows from the fact that the vertical axis is a 
crest and the horizontal one an anticrest. Inequality (R2d) is implied by the fact that the vertical 
axis is a valley and the horizontal one an antivalley. Inequality (R2e) is a consequence of the fact 
that the vertical axis is a crest and the horizontal one a valley. Inequality (R2f) fixes the orientation 
of the parameterization (cf. Section I4.1.4p . 
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